DARBOUX TRANSFORMATIONS OF JACOBI MATRICES 
AND FADE APPROXIMATION 



MAXIM DEREVYAGIN AND VLADIMIR DERKACH 



Abstract. Let J be a monic Jacobi matrix associated with the Cauchy trans- 
form F of a probability measure. We construct a pair of the lower and up- 
per triangular block matrices L and U such that J = LU and the matrix 
Zc = UL is a monic generalized Jacobi matrix associated with the function 
5c (-^) = ^F(X) + 1. It turns out that the ChristofTel transformation 3c of 
a bounded monic Jacobi matrix J can be unbounded. This phenomenon is 
shown to be related to the effect of accumulating at oo of the poles of the 
Fade approximants of the function "Sc although S'p is holomorphic at oo . 
The case of the (/L-factorization of J is considered as well. 



1. Introduction 

Let cr be a probability measure with an infinite support contained in K. Also, 
assume that all the moments of cr, i.e. 

t'da{t), j €1+, 

are finite. It is well known that the sequence {Pj}°^Q of monic polynomials orthog- 
onal with respect to a satisfies a three-term recurrence relation |lj 

(1.1) APj(A) = P,+i(A) + b,P,iX) + Cj_iP,_i(A), J e Z+ 

with the initial conditions 

P-i(A)=0, Po(A) = l, 

where bj e M and Cj > 0, j E Z+. Recall that the polynomials Pj are also called 
the polynomials of the first kind corresponding to cr. Besides the polynomials of 
the first kind, one can also associate with cr polynomials Qj of the second kind as 
solutions of (|l.ip with the following initial conditions 

Q_i(A) = -l, Qo(A) = 0. 

Clearly, the relation ()1.1|) can be rewritten as follows 

Jp{X) = Ap(A), 
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where p — (Pq, ^'i, ^2, • • • ) and J is a semi-infinite tridiagonal matrix of tlie form 

/bo 1 \ 

Co bi 1 

&2 



(1.2) 



J 



Cl 



V 



7 



The matrix J is said to be the monic Jacobi matrix associated with a. 

As is known [5], a monic Jacobi matrix J admits the unique i[/-factorization 
J — LU with the bidiagonal matrices L and U having the forms 



(1.3) 



L = 



/I 

h 



V 







u = 



1 

U2 



\ 



1 

M3 



if and only if the condition 
(1.4) 



p,(o)^o, jez+ 



is satisfied. T/ie Christoffel transformation of J is defined as follows 

J = LU ^ Jc ■■= UL. 

The matrix Jc is the tridiagonal matrix associated to the measure tda{t) |6]. 
On the other hand, if s_i is a real number such that 

(1.5) Q,(0)-s_iP,(0) ^0, jeZ+, 

then the monic Jacobi matrix J admits the C/L-factorization 

J^UL 

where U and L are defined by (|1.3p . The matrix 

Jg if/ 

is clearly a tridiagonal matrix and it is called the Geronimus transformation of 
J with the parameter s_i. Notice that the i[/-factorization of a monic Jacobi 
matrix is unique while the [/i-factorization depends on the free parameter s_i. 
Both the Christoffel and Geronimus transformation are also called the Darboux 
transformations [6] . The study of discrete Darboux transformations was originated 
in [25^ and has been further developed in [15], [29], [30]. Also, note that discrete 
Darboux transformations are applied in bispectral problems [TB] , Toda lattices ^2B] 
and numerical linear algebra [19]. We refer the reader to [6] for more details. 

In the present paper we extend the Darboux transformations beyond the con- 
ditions (|1.4|) and (|1.5|) . Namely, for an arbitrary monic Jacobi matrix J a pair of 
lower and upper triangular block matrices L and U is constructed such that J = LU 
and the matrix 

Zc^UL 

is a monic generalized Jacobi matrix in the sense of [llj . i.e. Zc is a tridiagonal 
block matrix with blocks specified below in Definition 12.41 
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Similarly, for any monic Jacobi matrix J and the free parameter s_i the UL- 
factorization J — il£, with the lower and upper triangular block matrices £ and il 
is found. In this case the Geronimus transform 

is also a tridiagonal block matrix (a monic generalized Jacobi matrix). 

Recall that the monic Jacobi matrix J can be also associated with the Nevanlinna 
function (see [1]) 

nx, ^ I pq. 

Analogously, one can associate the generalized Jacobi matrices 5c and with the 
functions and ^g, respectively (see [H]). We show that for the case in question 
they have the following forms 

;?c(A) - AF(A) + 1, ;Jg(A) = -^ + ^, 

where s_i is a free parameter. 

It turns out that the Christoffel transformation ^g of a bounded monic Jacobi 
matrix can be unbounded, that is, the entries of Zc are not necessarily bounded. In 
Section 6 we relate this phenomenon to the effect of accumulating at oo of the poles 
of the Pade approximants to ^c- In fact, we show that the L/7-factorization of J 
is bounded in the sense that L and U are bounded if and only if the diagonal Pade 
approximants to the corresponding function converge to 5c locally uniformly 
outside of a finite interval of the real line. A similar situation takes place in the 
case of the Geronimus transformation. 

The paper is organized as follows. In Section 2 we recall some basic facts and 
definitions. In particular, the definition of monic generalized Jacobi matrices is 
given. The Darboux transformations for arbitrary Jacobi matrices as well as for 
monic generalized Jacobi matrices are presented in Sections 3 and 4. Next section 
deals with triangular factorizations of arbitrary symmetric Jacobi matrices. In Sec- 
tion 6, we restrict our consideration to the case of probability measures supported 
on [—1, 1] and some criteria for the locally uniform convergence, outside of a finite 
interval, of the diagonal Pade approximants to the functions and ^g in terms of 
the LU- and [/L-factorizations are obtained. Finally, in Section 7 we provide the 
reader with some concrete examples. 

2. Preliminaries 

2.1. Classes D^^^ and DI^q. Let N be the class of all Nevanlinna functions which 
map the upper half plane C4. into the upper half plane C+. We will say (cf. [17]) 
that a Nevanlinna function F G N belongs to the class N_2n if for some numbers 
■Soi • • ■ J S2n G K the following asymptotic expansion holds true 

(2-1) ^(-^) = -T-S----Afe+''(A^)' 

where A^oo means that A tends to 00 nontangentially, i.e. inside the sector e < 
arg A < TT — e for some e > 0. Let us set 

n>0 
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Recall [1] that F e N_oo if and only if it admits the following representation 

FW = / ^, 

where tr is a bounded measure supported on the real line and Sj = j^Pda{t), 
j € It will be sometimes convenient to use the following notation 

oo 

(2.2) F(A)^-5]-^, A^oo 

to denote the validity of ()2.ip for all n G Z+ . The Jacobi matrix J associated with 
the probability measure a/sQ will be also called the Jacobi matrix associated with 

FeN_oo. 

Definition 2.1. Let us say that a function ^ meromorphic in C+ belongs to the 
class if it admits the asymptotic expansion 

(2-3) A^----' 

with some Sj G M (j G Z+) and \^{X) + Sq G N_oo- 

Definition 2.2. Let us say that a function ^ meromorphic in C+ belongs to the 
class D^Q^ if there exists a number s_i S M \ {0} such that 

In what follows we use the Gothic script for all the notations associated with the 
Di^j^-functions and the Roman script for the N_oo-functions to avoid confusion. 

2.2. The Schur transform of the D^^^-functions. Let be a nonrational func- 
tion from D^j^ having an asymptotic expansion of the form (|2.3p . Define the set 
Af{s) of normal indices of the sequence s = {Si}^o 

(2.4) = {n, : dct(s,+fe) ■ Uo ^0, j - 1, 2, . . . }. 

We will say that the function ^ g '^-oo is normalized if the first nontrivial 
coefficient in its asymptotic expansion ()2.3p has modulus 1, i.e. |Sni-i| = 1- 
The following statement can be found in [11] Theorem 4.8]. 

Theorem 2.3. Let ^ be a nonrational normalized function of the class D^^, let 
the sequence 5 = {Sj}°^Q be defined by the asymptotic expansion p.3p . and let 
A/'(s) = {rii}^]^ be a set of normal indices of s. Then ^ admits the expansion into 
the following P -fraction 

^ ' Po(A) - Pi(A) - ••• - PAr(A) ' 

where pt are monic polynomials of degree ii :— ni+i — rii (< 2), ct = ±1, b,; > 0, 
i € Z+, and Wq — 0. 

For the convenience of the reader we sketch the proof. Actually, the proof is 
based on the following step-by-step Schur process [11] (see also [1], [7]). If So 7^ 
then Ui = 1 and the polynomial po is defined as follows 

Po(A) = A - ^ 



So So 



So Si 

1 A 
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If 5o = then Si 7^ 0, ni = 1, and 
Po(A) = 



det(s,+j)i^-^o 

In both cases the Schur transform 3^ of e 

1 



Si S2 

Sl S2 S3 

1 A A2 



is defined by the equahty 



(2.6) 



= eoPo(A) + b^^(A), 



where eo — Sni-i and bo is chosen in such a way that ^ is normalized. Also, it is 
easy to see that ^ E D^oo (see [H] for details). Further, considering the asymptotic 
expansion of ^Ji := ^ 



one can construct the number ei — ±1, the monic polynomial pi, and the function 
i?2 := Continuing this procedure leads us to (|2.5p . 

2.3. Generalized Jacobi matrices associated with the D^j^-functions. The 
P-fraction p.Sp enables us to associate a generalized Jacobi matrix with the func- 
tion ^ € (see [8l[9l[Tg). 

Definition 2.4. Let be a nonrational normalized D^^^-function having the P- 
fraction expansion (j2.5p . Then a monic generalized Jacobi matrix associated with 
the function ^ e D^^^ is the tridiagonal block matrix 

/So Do \ 



(2.7) 



£1 »2 



where the diagonal entries have the following form 



Po ' 



Po 



if«,= 

if h = 



1; 

2, 



the super-diagonal is built up with ij x tj+i matrices 

^0^ 



and the sub-diagonal consists of tj+i x tj matrices 



here Cj — ejCj+ib^. 

Remark 2.5. Actually, one can associate a monic generalized Jacobi matrix Z with 
an arbitrary nonrational function 5^ € Di^,^. Namely, every function ^ e D^,^ can 
be normalized as follows ^nor = i?/|Sni-i| and the monic generalized Jacobi matrix 
associated with Jnor will be also called associated with ^ (see Remark 2.5). 
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Setting C-1 = eo, let us define polynomials of the first kind ^j(A), j £ Z+, as 
solutions Uj — *Pj(A) of the following system: 

(2.8) <^j-iyj-i - PjWvj + Vj+i = 0, j 

with the initial conditions 

(2.9) 2/_i =0, yo = 1. 

Similarly, the polynomials of the second kind 0j(A), j € Z+, are defined as solutions 
Uj = nj(A) of the system ()2.8p with the initial conditions 

(2.10) y-i^-1, yo = 0. 

It follows from (12. Sp that is a monic polynomial of degree rij and Qj is a 
polynomial of degree xij — £o with the leading coefficient eo- The equations ()2.8p 
coincide with the three-term recurrence relations associated with P-fractions |24| 
(see also [18l Section 5.2]). The following statement is immediate from (12. 8[) . 

Proposition 2.6. f[10j). Polynomials and (Xlj and Qj^i) are coprime. 

Recall (see [21]) that the polynomials can be found by the formulas 

/ So Si ... Snj \ 



(2.11) <Pj(A) = — det 



and hence 
(2.12) 



Srij - 1 

V 1 



*P.-(0) - ^^det(s,+,+i)"Co' ^^(0) = 



det 



/ 

So 



det(si+fe)"^fc^o), 



So 

Si 



Srij-l 
Sn, 



S2nj-1 / 



since Oj(A) = "^t ( 



%(A)-q;,(t) 



, where £f is a linear functional such that 



Introducing for arbitrary j G Z_|_ the shortened matrices 

/So So \ /»! S)l 



(2.13) 3[o,,] = 



£0 «i 



V 



£1 *B2 



V 



one can obtain the following connection between the polynomials of the first and 
second kind Qj and the shortened Jacobi matrices Z[o,j] and (for the 

classical case see [5] Section 7.1.2]). 

(2.14) q3j(A) =det(A-a[oj„i]), Oj(A) = eo det(A - ^[ij-i]). 
Remark 2.7. Let us define an infinite matrix G by the equality 

e„ iftj^l; 

(2.15) G = diag(Go,...,G„,...), G, 



if e, = 2, 
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and let i'^^ (G) be the space of i?^-vectors with the mner product 



(2.16) 



[x,y] = {Gx,y)g 



x,y e£l 



[0,00)' [0,00)- 



The inner product (|2.16p is indefinite, if either kj > 1 for some j e Z+, or at least 
one Ej is equal to —1. The space oo)(^) equivalent to a Krem space (see [3]) 



if both G and G are bounded in 



[0,oo) 



If k. 



1 for all j big enough, then 



^fo oo)(^) ^ Pontryagin space. 

The m-function of the shortened matrix Z[a.j~i] is defined by 
(2-17) m[oj-i](A) = [iZio,j-i] - A)-ieo,eo], 

where cq = (1, 0, 0, . . . , 0)^ G C"^ . Due to (|2A4)) it is calculated by 

det(A _ Qj{X) 



(2.18) 



m 



-i](A) 



det(A-3[oj_i]) *P,(A)- 



Remark 2.8. Let us emphasize that the polynomials *Pj and Qj have no common 
zeros (see [31 Proposition 2.7]) and due to (12.181) the set of holomorphy of m^oj-i] 
coincides with the resolvent set of ^[o.j-il- 
It was shown in [9 that m[oj-i] has the following property 

(2.19) ^(A)-m[o,,_i](A)=o(^^^5i^), A^c^. 

The latter means that one can reconstruct the sequence {Sfej^g by the monic 
generalized Jacobi matrix Z- Namely, it follows from (|2.19p that for any k € Z4. 



(2.20) 



Sfe 



(^[oj-i]) 



eo,eo 



for all j such that n,- > 



3. The Christoffel transformation and its inverse 



From now on we always assume that F G N_oo (or ^ G D^^^) is a nonrational 
function. Then there is a monic Jacobi (or generalized Jacobi) matrix associated 
with F (or 5^, respectively). 

3.1. L[/-factorizations of Jacobi matrices. In the following proposition we in- 
troduce a block L[/-factorization of an arbitrary Jacobi matrix. 

Proposition 3.1. Let J be a monic Jacobi matrix associated with F G N_oo with 
the asymptotic expansion (|2.ip . let {nj}°°^^ be the set of normal indices of the 
sequence s — {Sj}°°^Q :— {sj+i}j^o and let tj :— rij+i — rij, j G Z+, where rig = 0. 
Then J admits the following factorization 

(3.1) J = LU, 

where L and U are block lower and upper triangular matrices having the forms 
(h„ \ (Uo Do \ 



(3.2) 



L = 



L2 







\ 



U = 



■J 







U2 



■J 
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in which the sub-diagonal of L consists ofij+i x matrices 



(3.3) 



L 



the diagonal entries of U are of the form 



(3.4) 



1 



and t/ie super- diagonal of U is built up with tj x ij+i matrices 

(3-5) D, = (J (J), (1 0), (1). 

Moreover, the factorization (j3.ip is unique and the following relations hold true 



(3.6) 



,0-) 



6n-+i, j e 



P„,(0) ' ' 

Proof. It is natural to consider the following relation 

/ Uo Do \ /bo 1 



LU 



L2U1 L2D0 + U2 



V 



Co 61 1 

ci 62 



V 



which enables us to find the entries of L and U. By the definition, there are two 
possibilities for the matrix Uo- If to = 1 then 

A(0) PnAO) 



bo 



Po{0) Pno(O)' 



In case to = 2 we have 







1 



(0) „,(o) 







60 1 
Co 61 



and, hence, one obtains the equalities 



(0) (0) , 

u'o = Co, = bi. 

Since to = 2 then si = 0, Pi(0) = and (HH) yields 



,(0) 



Co 



^2(0) P„,(0) 



PoiO) Pn„{0) 

The equality 60 = follows from the relation bo — ^ (see [1]). 
Now, assume that 

We will consider 4 cases. 

Case 1. tj — tj-i = 2, j e N. In this case one obtains the following relations 
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which immediately imply 



(3.8) 



.0) 



1 

&n,+l 



It follows from = = 2 that det{s,+i)"]tl = 0. Hence by ([^1^ Pn,+i{0) 
-Pnj-i(O) = and then yields 6„^. = 0. Moreover, from (ll.ip one gets also 



(3.9) 



Pn,+2(A) + (5„^.+l-A)P„^. + l(A) 



^n,(A) 



A=0 



Pn, + .(0) 
^'n,(0) ■ 



Case 2. L = 2, L_i = 1, e N. In this case we have 





f) 


h 


1 


4^) 









6n, 1 
Cn, &n,+l 



Since Pnj+i(0) = 0, one concludes as in 

(3.10) u'^l^ = bn^+, 

In order to see that the relations 
(3.11) 



that 



0) ^n.+i(0) 



hold true simultaneously, let us set Ij :— bn ■ Thus, according to the assump- 
tion p.7p . the second relation in p. lip can be rewritten as follows 

(3.12) &n,Pn,(0)+C„^-lFn,-l(0) =0, 

and formula p.l2p follows from ()l.ip and the fact that Pn^+i(0) = 0. 

Case 3. ij = 1, {^-i = 2, j e N. The identities Ij = 1 and Ij^i = 2 yield 

i,[/,_i = (0 /,)-(0 c„^_i), L,D,^, + U,=u^i^ =h^,^. 

Consequently, by taking into account the equalities Pnj -i(O) — and n^+i = + 1 
we get that 

(3.13) 



^n, + ,(0) 
^n,(0) ■ 



Case 4- ij = I, ij-i — 1, j E N. This case reduces to the following relations 



(3.14) i,C/,_i = lA''' = cn,-i, LjD,., + U, = + 



Eliminating from the first equation in (|3.14p and substituting it to the second 
one, we get 

Pn, (0) + Cn^- - 1 Pn, - 1 (0) _ P„,^, (0) 



(3.15) 



,(]) 



^n,(0) 



^n,(0) 



□ 



Theorem 3.2. Let J he a monic Jacobi matrix associated with F € N_oo and let 
J = LU be its LU factorization. Then the matrix ^ — UL is the monic generalized 
Jacobi matrix associated with XF{X) + sq G Dl^. 
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Proof. The fact that J = UL is a monic generahzed Jacobi matrix can be easily 
verified by straightforward calculations. 

We can assume that the function F is normalized and let F have the asymptotic 
expansion (|2.2p . Then the function XF{X) + 1 has the asymptotic expansion 

AF(A) + l^-^-g-... (A^oo) 

and is not necessarily normalized. Since according to (|3.6|) we have 

(o)__ PnM-i ^osi, if«o = l; _ r if«o = l; 

°° - 'V„„(0) "1 eoS2, if«o = 2 -\\S2\, if«o = 2 

the function ^nor = iF{X) + l)/eoUg'^'' is normalized. 
Further, observe that 

(3-16) ^ioj]^Q = ^0, Glo.j]eo = -%;^[o.j]^o, j e 

where the shortened matrices i[o.j] , ^[Oj] > ^-nd G[oj] are defined analogously to (|2.13l) . 
Clearly, for j big enough 

where J[o,n is the leading principal submatrix of J constructed from the first 

rij+i rows and columns of J. Thus, we have 

Sn = (eo,^[o,n^+i_i]eo) = (^ep, £[oj][/[oj] . . . Ljo^jjC/jo^jj^ eoj = 

n times 



n times 

n e 

n— 1 times 

Further, using f|3.16p gives the equahty 

Sn = eo4°^ (eo, ([/[oj]^[Oj])"~^G'eo) . 
Hence for sufficiently large j one calculates the moments 

Sn 



(3.17) 



{dJo,j]) eo,eo 



(0) ' 



which coincide with the coefficients in the asymptotic expansion of Snor- In view 
of (|2.20p this means that the generalized Jacobi matrix ^ is associated with the 
function Ai^(A) + sq e Dl^. □ 

The transform J = LU t-^ Z = UL is called the Christoffel transform of the 
Jacobi matrix J. 

Corollary 3.3. Let J be a monic Jacobi matrix associated with F G N_oo with 
the asymptotic expansion (|2.ip Assume that Pj{0) ^ for all j G N. Then Wj ~ j 
for J G N and 

(i) The matrix J admits the LU -factorization (|3.ip . where L and U are lower 
and upper triangular matrices having the forms p.2p with 

(3.18) D,=l, [/^- = _i2±iM, = - (7,+!, JGZ+. 
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(ii) The matrix ^ = UL is the monic tridiagonal generalized Jacobi matrix 
associated with XF(X) + sq G DI^q. 

The statement of this corollary is contained in [6], where, in fact, more general 
monic three-diagonal matrices of the form ()1.2|) were considered (their matrices 
are associated with a class of quasi-definite linear functionals which includes finite 
measures on the real line as a subclass). 

3.2. [/L-factorizations of generalized Jacobi matrices. Here we present a 
block C/L-factorization of monic generalized Jacobi matrices. 

Proposition 3.4. Let ^ be a monic generalized Jacobi matrix associated with ^ G 
T^Zoo! 1st {j^jl^i be the set of normal indices of the sequence s — {Sj}°°^Q defined 
by the asymptotic expansion (|2.3p of the function ^, tj :— rij+i — rij, j G 
where Uq — 0, let s_i G M 6e the same as in Definition \2.2l let ^j(A) and Qj{X) 
be polynomials of the first and the second kind associated with the sequence 5 — 
{5j}j°^Q, and let the polynomials *Pj be defined by 

$,(A)==*p,(A)--^Q,(A). 



S-l 



Then: 



(i) ^,{0)^0 for J eZ+; 

(ii) the GJM 3 admits the following factorization 



(3.19) 



(3.20) 



where U and L are block lower and upper triangular matrices having the 
form (j3.2p - (13.51) . Moreover, the following relations hold true 



,(J) 



lj+1 - 



^ , ilrx 

*P,(0) 



1 

S-l ' 



ij+i 



Proof, (i) Notice first that the sequence {sj}°^Q — {Sj-i}°^Q satisfies the asymp- 



totic expansion 



of the function F{X) = 



det 



A-i 

So 



So 
Si 



A 

Sk \ 



G N_ 



Then 



S/c+l 



^0, fc G 



V S/c Sfe+i ... S2k+l/ 

by Definition 12.21 and the property of Nevanlinna functions. In particular, one gets 
from (|2.12p the relation 

f 

So 



Sl 



Sn^-i =s_idet 



det 



S2nj-1 y 



So 
Sl 



Srlj-l 
S„, 



Vsn,-1 



S2n,-1/ 



- i-iy+' (s-i^,(0) - 0,(0)) detis.+k)Z=l + 0, 

which shows that (0) ^ for j G Z+ . 
(ii) Further, the equality Z — UL yields 



(3.21) 



(i) 



,0') 



^3' 
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Setting Mg°' := l/s_i, we see that 



ti - -Po - 



1 



^i(O) 



s-i *Po(0) 

Next, the second formula in p.20p follows by induction with the help of the 
second and third equations in p.2ip . and the fact that the polynomials sat- 
isfy (EH). □ 

Remark 3.5. It should be noted that, actually, both the matrices in the J7L-decom- 
position p.l9|) depend on s_i, that is, U = C/(s_i) and L = L(s_i). 

Theorem 3.6. Let ^ be a monic generalized Jacobi matrix associated with ^ G 
T)Zoo, Is-t s_i G R &e the same as in Definition \2.S[ and let ^ ^ UL be its UL 
factorization of the form p.l9p . p.2|) - p.5p . Then the matrix J = LU is the monic 
Jacobi matrix associated with — ^-j^ G N_oo- 

Proof. The fact that J is a classical Jacobi matrix can be easily verified by straight- 
forward calculations. The rest of the proof can be done by reversing the reasoning 
given in the proof of Theorem 13.21 □ 



4. The Geronimus transformation and its inverse 
4.1. LJ7- factorizations of generalized Jacobi matrices. 

Proposition 4.1. Let ^ be a monic generalized Jacobi matrix associated with ^ G 
D^oo '^'^'^ '■— ^j+i ^'^j ■ j ^ ^+-' where tIq = and {nj}°°^i is the set of normal 
indices of the sequence s — {Sj}°^o defined by (|2.3|) and let ^j(A) be polynomials 
of the first kind associated with the sequence s = {Sj}°°^Q. Then *Pj(0) 7^ for all 
j G and the GJM Z admits the following factorization 

(4.1) a = £11, 

where £ and it are block lower and upper triangular matrices having the forms 



(4.2) 



/(So 
£1 



V 





Ci 





^2 



u 



/ito Do 
ill Di 

iio 



\ 



in which the sub-diagonal of £, consists of ij x t 

0\ fo 

[, r I ( 



(4.3) 



j-i matrices 



(0 



■3/ VJ/ 

1)j are of the form (j3.5p and ilj , are tj x tj matrices 



,(j) 



(4.4) 



ifh = 2. 



1, 



Moreover, the following relations hold true 



(4.5) 



< — ■> 



1 



Pi ' 



= 1; 
=2. 
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Proof. First, notice that it follows from formula ()2.1ip that 

det(Si+fc+i)"i_n 

since ^ e D^'^g^ and Si := s^+i {i E Z_|_) are moments of the Nevanlinna function 
Consider the product £il of the matrices £ and il 



£iilo £i2)o + £iUi 

£2iii £22)1 + €2112 



V 



/ V 



/'Bo 


So 




\ 






©1 






£1 






V 






•7 



Comparing it with the matrix ^ in (12.71) one finds the entries of £ and il. 
Next, we will prove formula (|4.5I) . 
If £0 = 1 then the equality €oiio = ®o yields 



(4.6) 

In the case = 2 we have 
and, again, we obtain the equalities 











eo J 


-{ 



Po 



Now, assume that 
(4.7) 



(0) _ *Pi(0) 
*P.(0) 



-pr- 



U,_i 



^.-i(0) 



We will analyze four cases. 

Case 1. Ij = — 2, j G N. In this case one gets the following relations 



(4.8) 




[jUj_i 





c,_i or 



(4.9) £,D,_i + 

Hence one obtains 



1 



Po 
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Case 2. Let ij = 2, tj-i = 1, j G N. Then gS]) holds true and (|48j) takes the 
form 

■ \ / 

Hence one obtains (|4.10p and (|4.11l) . 

Case 3. Let «j = 1, tj^i = 2, j e N. Then 

(4.12) = ([,u,_i 0) = (c,_i 0), 

(4.13) £,D,_i + = [, + u, = -p(^') . 

Case 4- In the case tj — 1, tj-i = 1, j G N the equahty (j4.13p holds true 
and (|4.12p takes the form 

'j"j-i = S-i- 

In both cases the calculations in (|4.1ip are still in force. □ 

Theorem 4.2. Let ^ be a monic generalized Jacobi matrix associated with ^ S 
Dj^g^ and let J = £il be its LU -factorization of the form (j4.ip - (j4.5p . Then the 
matrix J ~ 5I£ is the monic Jacobi matrix associated with A5'(A) + Sq G N_oo- 

Proof. The fact that J — U£ is a monic Jacobi matrix follows by straightforward 
calculations. Next, notice that 

(4.14) ^io,]]^o = eo, ilfojiGfojieo = aeo, j e1+, 
where 

eoUo, if?o = l; 



a 



eo, if«o-2 



and the shortened matrices £[o j] , ^[o,j] i ^nd G[o j] are defined analogously to (j2.13p . 
We can assume, without loss of generality, that the function ^ is normalized. Then 
it follows from (|4.6p that a = Si both for to = 1 and {q = 2. 
It follows from (|2.20p and (|4.ip that for j big enough one has 

Sk = (eo,5fo,j-i]G[oj-i]eo) = (ep, £[o,j]ii[oj] . . . £.[o,j]H[o,j] G[o^j]eo 

k times 

£^eo,il[oj]-C[oj] ■ ■ •^[Oj]-C[oj]iI[Oj]G[oj]eo), k eN 



k—l times 

Further, using (I4.14p we get 

Sk = a (eo, (ii[oj]£[o,j])''"^eo) , 
which for sufficiently large j can be rewritten as follows 

This implies that the Jacobi matrix J is associated with the normalized function 
F(A) = Mm±£o e N_oo and, thus, also with AS-(A) + sq. □ 
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4.2. t/i-factorizations of Jacobi matrices. 



Proposition 4.3. Let J be a monic Jacobi matrix associated with F G N_oo which 
has the asymptotic expansion (j2.1|) . s_i G M, and let Pj{X) and Qj{X) be polyno- 
mials of the first and the second kind, respectively, associated with the sequence 
s = {sj}°^o. Then: 

(i) the normal indices rij {j G N) of the sequence s = {sj-i}j^o '^'^'^ char- 
acterized by the conditions 

i^n,-i(o)^o, jeN, 

where Pj{X} := Qj{X) - s_iP,(A); 

(ii) the Jacobi matrix J admits the following factorization 



(4.15) 



(4.16) 



J = it£, 

where £ and ii are block lower and upper triangular matrices having the 
form (|4.2p - (j4.4p . Moreover, the following relations hold true 



U 



^n, + .-l(0) 
^^n,-l(0) 



Proof, (i) The normal indices xij {j e N) of the sequence s = {sj_i}°^Q can be 
characterized by the conditions 

/s-i 



dj :— det 



So 



So 

Sl 



Sn,-2\ 
Sn,-1 



7^ 0, j e N. 



\Snj-2 Suj-l ■■■ S2nj-3/ 

Now the first statement follows by the equalities 

/ So ... Snj.-2^ 

=,s_idet' ' ' '° '"'-1 




\Sn,_l Sn,+1 

= (-1)"^ (s_iP„^_i(0) - Qn,-i(0)) det{s,+k)::,Zl + 0, 
(ii) Comparing the matrix 

/itoeo+So£i So^i \ 



S2nj-3/ 



il£: 



ill£l 



ili£i+2)i£2 2)i£2 



with the matrix -3 in (|1.2p one finds the entries of £ and H. Let us consider four 
cases. ^ ^ 

Case 1. to = h = 2. In this case s_i = Po(0) 0, Pi(A) = 1 and P2(A) 
02(A) = A- fei.^ 

Since hi = —P2{0) =0 the equation 



(4.17) 



ilo£o + 2)o£i 



eo 1 
uo 



&o 1 

Ci 61 
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is solvable and it follows from 

^0 h\ fO C2 



(4.18) lli£i = 

that 

h = C2 = - 



0/ \0 

^3(0) 



A(0) 

Case 2. Iq = 2, £i = 1. In this case s_i = Po(0) = 0. It follows from the 
equation 



V"o hj Oi 
that 

Case 5. In the case to = 1, Ji = 2 one gets Po(0) 7^ 0, Pi(0) = 0, and the 
equations take the form 

iloeo +2)0-21 =Uo = 60, 



Hence 



[ 



, ^^2(0) 
II = Ci = - ^ . 

Po{0) 

Case 4- In the case = 5i = 1 one has s_i — Po{0) 7^ 0, and the equation 
Uo£o+2)o£i =Uo + [i =60, 
contains a free parameter Uq. Setting Uq = one obtains 

[1 = 5o - uo = • 

Assume now that ()4.16p is satisfied for some j G N and consider again four cases. 
1) ij — ij^i = 2. These conditions can be rewritten as follows 

Pn,.l(0)7^0, P„^. + i(0)7^0, P„,(0)=0, Pn^.+2(0) = 0. 

This implies bnj+i = and hence the equations 

(4.19) it,., + ^^J, 



_ fO c„^.+i 



^1^1 -{0 "0 J-[o 

are solvable and 

(4.20) = = -W^^^ 

Pnj + i-l(U) 

2) tj^i = 2, tj = 1. In this case 

Pn,_l(0)^0, P„,(0)=0, Pn^+i(0)^0, P„^.+2(0)7^0. 



DARBOUX TRANSFORMATIONS AND FADE APPROXIMATION 



17 



It follows from the equation 
that 

3) In the case ij-i = 1, tj = 2 one obtains 

Pn,_l(0)^0, Pn,(0)^0, P„^ + i(0) = 0, P„^.+2(0)^0. 

It follows that 
Hence 

i^-^^j Ij + l - Cn,+1 - — -. 

Pnj + i-l(,Uj 

4) In the case = ij = 1 one has 

P„^._i(0)^0, P„,(0)^0, P„^+i(0)^0, 

and the equations 

(4.23) %(£^- + Dj£j+i = Uj- + [j+i - 6n, , 

yield 

Cn,-1 



(4.24) 



&„,P„, (0) + C„^.-lP„,-l(0) _ Pn, + .-l(0) 



J^n,(0) Pn, + ,-l(0) 

This completes the proof of ()4.16p . □ 

Remark 4.4. It should be noted that, actually, both the matrices in the t/L-decom- 
position (|4.15p depend on s_i, that is, il — H(s_i) and £ = £(s_i). 

Theorem 4.5. Let J he a monic Jacobi matrix associated with F G N_oo, S-i G 
and let J ~ il£ be the corresponding U L factorization of J of the form (j4.2p - (|4.4I) . 
Then 2 — ■Gil is the monic generalized Jacobi matrix associated with ^-^^ — G 

Proof. The fact that Z is the monic generalized Jacobi matrix can be easily verified 
by straightforward calculations. The rest of the proof can be done by reversing the 
reasoning given in the proof of Theorem 14.21 □ 



The transform J = U£ J = £il in Theorem 14.51 is called the Geronimus 
transform of the Jacobi matrix J with the parameter s_i. 
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4.3. Jacobi matrices associated with generalized Stieltjes functions. Let 

K be a nonnegative integer. Remind that a function F, meromorphic in C \ R, is 
said to belong to the class if the domain of holomorphy p{F) of the function F 
is symmetric with respect to R, F{X) = F{X) for A G p{F), and the kernel 

\Nf{X,X)^F'{X), XepiF) 

has K negative squares on p{F). The last statement means that for every n € N 
and Ai, A2, ■ . ■ , A„ G p{F), n x n matrix {^pi^i, j=i has at most k negative 
eigenvalues (with account of multiplicities) and for some choice of n, Ai, A2, . . . , A„ 
it has exactly k negative eigenvalues (see [31] )■ Clearly, No = N. 

Definition 4.6. Let us say that a function F holomorphic in C+ belongs to the 
generahzed Stieltjes class S^'' if F{X) e N and X^^F{X) e N^. Let us set 



—00 



According to [T31 Theorem 2.4] any function F E S'l^ admits the integral 
representation 



^ A- /■= 



da{t) 



(4.25) 

where Aj > (j = 1, . . . ,p), p = k and cr is a finite measure on [0, +00), such that 

/•OO 

(4.26) / f-^dait) < 00 for ah n € N. 

Jo 

Similarly, any function F € SZ'^ admits the integral representation (|4.25p . where 
Aj > (j ~ 1, . . . ,p), a is a finite measure on [0, +00), which satisfies (|4.26p and 



(4.27) p : 



if < i^(O-) < 00; 
K, ifF(0-)<0. 



This implies, in particular, that every function F E belongs either to or 
to 8"+!. 



Corollary 4.7. Let F E have the asymptotic expansion (j2.ip . let J be a 

monic Jacobi matrix associated with F , let {nj}"^^ ''^^ 0/ normal indices 

of the sequence s — {sj+i}°^q, and let J = LU he its LU -factorization of the 
form ^Jl-I^. Then: 

(i) the sequence tj := rij^i — xij, j E is stabilized, i.e. there is N > such 
that tj = 1 for j > N; 

(ii) the matrix J = UL is the monic generalized Jacobi matrix of the form 



^i21 i)[Af+l,oo 

■II[Ar+i,oo) a monic Jacobi matrix, -3[o,Af] a monic generalized Jacobi 
matrix of the form (j2.7p ; 
(iii) the matrix J = UL is associated with XF{X) + sq E Dl^. 
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5. Generalized Cholesky decomposition 

Recall that the Cholesky decomposition is a decomposition of a symmetric 
positive-definite real matrix A as the product of a lower triangular matrix L and 
its transpose L^, that is A = LL^ . The decomposition can be also rewritten in 
the form A — LAL^ , where A is a diagonal matrix and all the diagonal entries of L 
are equal to 1 . In this section, we give a generalization of the latter decomposition 
to an arbitrary symmetric Jacobi matrix J^. 

At first, note that every monic Jacobi matrix J can be reduced to a symmetric 
Jacobi matrix via the following transformation 

(5.1) Js = 

where 5" — diag (1, y/co, ^/coCl, ^cqC\Ci, . . . ); here the square root is chosen to be 
positive. In fact, every symmetric Jacobi matrix can be represented in such a way. 
So, now we can reformulate Proposition 13. II in an appropriate way. 

Proposition 5.1. Let Js he a symmetric Jacobi matrix associated with F G N_oo 
and let tj := rij+i — rij, j G Z,-^, where no = and is the set of normal 

indices of the sequence s ~ {-Sjlj^o defined by the asymptotic expansion (|2.3p of 
the function AF(A) + sq. Then J admits the following generalized Cholesky decom- 
position 

(5.2) Js=LAL^, 

where L is a block lower triangular matrix having the form 

\ 



in which the sub-diagonal of L consists of ij+i x tj matrices 

2.+i = (V o)' (V)' (^"^^^ 

and A — diag (Aq, Ai, . . . ) is a block diagonal matrix with the entries 

Moreover, the factorization (j5.2p is unique. 
Proof. Making use of (|3.ip and (|5.ip gives 



L = 



A, 




(5.3) 



Js ^ ^'"'Lf/*. 
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Further, let = diag (v^o, 5*1, ^2, . . . ) be a partition corresponding to the one of 
L. Then (15.31) can be rewritten as fohows 



(5.4) 



Co 



v 



Co 
61 



/Cl 

b2 



/A 

ZiAo 



LiBo + Ai 
L2A1 



L2D0 + A2 



•7 



where 

Since the matrices in 



are symmetric we have that Aj — Kj and Dj 
(Lj+iAj)^ = KjLj^^ for j e Z+. This observation and the representation ([57 
yield (jOl). □ 



6. Convergence of Fade approximants for D^^ 



FUNCTIONS 



In this section we restrict our consideration to the case of probability measures 
supported on [—1, 1]. As a consequence, the corresponding monic classical Jacobi 
matrices will be bounded. To specify the results of this section, let us recall the 

notion of Fade approximants. Suppose we are given a formal power series F{\) = 

00 

' with Sj E M. Let L, M be positive integers. Then an [L/M] Fade 



E 



approximant for F is defined as a ratio 



4 [L/M] (1 



of polynomials of formal degree L and M, respectively, such that 

5[L/A/](o) ^ and 

L+A/-1 



(6.1) 



E 



+ J7[VM](A) = 



1 



A 



00. 



More information about Fade approximants can be found in 4 . In what follows C 
is considered as C U {00} equipped with the spherical metric. 

6.1. Pade approximants for Dl^^ functions. In a previous work we proved the 
following result. 

Theorem 6.1 (111^). Let a he a finite nonnegative measure on E = [—1, a] U [/3, 1], 
G [a, /?], and let 

tda{t) 



(6.2) 



t-X 



Then ^ G ^^00 '^^'^ sequence {ijf"^^"^i}°^Q converges to ^ locally uniformly in 
C \ ([~1 ~ Of] U [/?, 1 + e]) for some s > if and only if the condition 
lao\ -Pnj+l(0) 

(6.3) sup —;i — < 00 



is fulfilled for the polynomials Pj orthogonal with respect to a. 
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Remark 6.2. In [11, Section 5.3] we gave some sufficient conditions for the measure 
(7 to possess the property ()6.3|) . 



This resuh can be interpreted in terms of the Darboux transformations as fohows. 



Theorem 6.3. Let ^ G ^-oo have the form (|6.2p . Then the following statements 
are equivalent: 

(i) for the number 

S_i := / d(j{t), 

J E 

the operator U = C/(s_i) in the factorization (|3.19|) of Z corresponding to 
g' is hounded in 1"^^ ^y, 

(ii) the operator U in the factorization p.ip of J corresponding to the function 



F{X) 



da{t) 
t-X 



is bounded in i'^^ 



(iii) the sequence {^^"'^ ^"'^}'j^Q of diagonal Fade approximants converges to ^ 
locally uniformly in <C \ ([—1 — e, a] U [/?, 1 + e]) for some e > 0. 

Proof. Clearly, the first statement is equivalent to the second one by construction. 
So, it remains to prove the equivalence of (ii) and (iii). First, suppose that (iii) holds 
true and, therefore, from Theorem 16. II we see that ()6.3|) is fulfilled. Then, according 
to (|3.8p . p.lO|) . (|3.13p . and p.l5|) . the operator U is bounded. The implication from 
(ii) to (iii) follows from formulas p.l3p . p.lSp and Theorem 16. II □ 

Remark 6.4. Under the condition of Theorem lO from (|3S]), (|3TT0l) . ([3lT|) . (|3?T3l) . 
and p.l4|) we see that the operators U and L, given by (|3.1|) . are bounded or 
unbounded simultaneously. 



In some cases, we can specialize Theorem [ 
Corollary 6.5. Let ^ admit the following representation 

where a is a finite nonnegative measure on [0, 1], tj < 0, and a.j > for j — 1, . . . ,p. 
Then the following statements hold true: 

(i) the sequence {5^'"^^"^i}^o '^f diagonal Fade approximants converges to ^ 
locally uniformly in C \ ([0, 1] U {tk}l.^i); 

(ii) the polynomials Fj orthogonal with respect to a + X]j=i %©(■ ~ tj) possess 
the following property 

(here Q denotes the Heaviside step function); 

(iii) for the number 



sup 



< oo 



1 P 
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the operator U ~ U{s-i) in the factorization p.l9p of ^ corresponding to 
g' is bounded in £|j 

(iv) the operator U in the factorization p.ip of J corresponding to the function 
is bounded in i'^^ . 

Proof. Clearly, ^ E ^^oo- ia.ct, statement (i) for the class in question was 
proved in [27]. Statement (ii) is an immediate consequence of statement (i) and 
Theorem 16. II The rest follows from Theorem 16.31 □ 



Remark 6.6. As was mentioned in Subsection 4.3 the function F in (j6.5p belongs 
to the class S~p. 

6.2. Pade approximants for D^oo functions. In this section we present con- 
vergence results for diagonal Pade approximants to D^^^ functions. 

By re-examining 10, Theorem 4.16], [11] Theorem 3.3] and [TT] Theorem 5.5], 
we arrive at the following more general result. 

Proposition 6.7. Let 

da{t) 



F{X) 



1 i - A' 



Then for any sequence {tj }^]^ of real numbers the following relation holds true 



3 = 



(6 6) F["./n.] ._ Qn,(A)+r,Q„^._i(A) _ ^ 
^ ■ ^ ■ P„^.(A)+r,P„^._i(A) ^ A^+i ■ " VA^"^ 

where Pj and Qj are polynomials of the first and second kind, respectively, corre- 
sponding to a. Moreover, the sequence of the modified Pade approximants pl'^j/^jl 
converges to F locally uniformly in C \ ([— 1 — £, 1 -|- e]) for some e > if and only 
if the condition 

(6.7) sup |tj| < oo 

jeN 

is satisfied. 

Proof. The relation (16.61) is an appropriate reformulation of the results of [T] Sec- 
tion 1.4.2]. Next, notice that the symmetric Jacobi matrix J' associated with tr is a 
bounded linear operator in l'^^ Furthermore, using the classical versions of for- 
mulas (I2.14p . (|2.17p . and (|2.18p . the modified Pade approximant can be rewritten 
as follows 

i?[n./n.l(A)^ Q„,(A)+r,g„^_i(A) _ det(A- j(;;)_^,) 



Pn, (A) + T,P„^. _i (A) dct(A - ) 
= ((^[K-i]-A)-Ve), 
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(t) 

where JTjq is a rank-one perturbation of i7[o,nj-i] of the form 



[0,n,-l] 



V 



= J[omj-i] +diag{0, ...,0,Tj} 



B 2 bn^~l+TjJ 



(cf. [28l Proposition 5.8]). Now, let us suppose that (|6.7p is satisfied. Then one 
obtains 

(6.9) < ll^[o,n,-i]ll + |t,| < 1 + e, J e Z+, 

for some e > 0. It foUows from (|6.8I) . the inequahties 



1 



and (ESI) that 



for |A| > 1 + e. The pointwise convergence f["j/''^1(A) to F(A) for A G C\R follows 
from the proof of [281 Proposition 5.3]. Finally, it remains to apply the Vitali 
theorem. 

Let us prove the necessity by proving the contrary statement. So, suppose that 



sup I r,- 1 

jeN 



oo. 



Further, note that the poles of the modified Pade approximant coincide 



with eigenvalues of the matrix >7[Q^n^. Taking into account that i/[ori^ -il 
self-adjoint matrix, one obtains 



[0,n,~l] 



|AmQa;(J'[o n^^_l])| 



[0,nj-l]l 



^ l('7'[fn,_lien,^l,en,-l)| = |&n,-l +T, 



where \max{J\^Q[^._i-^) is the eigenvalue of the matrix J\^^'^._it^ with the largest 
absolute value. Since the sequence {&fc}^Q is bounded, we have that infinity is 
an accumulation point of the set of all poles of the modified Pade approximants 

Theorem 6.8. Let a he a finite nonnegative measure on [—1,1] and let 
(6.10) I 



1 r da{t) s_i 



where s_i is a fixed real number. Then ^ G ^-oo '^^d, the sequence {d^"'^^"'^^}j°^Q 
converges to ^ locally uniformly m C \ ([— 1 — £, 1 + e]) for some s > if and only 
if the condition 

Qn,(0)--s_iP„^.(0) 
(6.11) sup — -TTT. 7^ — — < oo 



sup 



)„,_i(0)-s_iP„^._i(0) 
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is fulfilled for the polynomials Pj and Qj of the first and second kind corresponding 
to a . 

Proof. Observe that by setting 

QnAO) - S-lPnAO) 



.i(0)-s_iP„._i(0) 



we have that 5^1"j/"^1(0) = s_i. Thus, due to (|6.ip and (|6.6p one has the followmg 
representation of the diagonal Pade approximants to 5^ 

Now the statement is an immediate consequence of Proposition 16.71 □ 

The above resuh can be reformulated in terms of the Darboux transformations 
as follows. 

Theorem 6.9. Let ^ G Dl^^ have the form (j6.2[) . Then the following statements 
are equivalent: 

(i) the operator £ in the factorization (j4.ip of J corresponding to 5" is bounded 

(ii) the operator £ = £(s_i) in the factorization (14.151) of J with the parameter 
s_i = So corresponding to the function 



da{t) 



X 



'[0,oo)' 

(iii) the sequence {^^"'^^"^^}j°^Q of diagonal Pade approximants converges to ^ 
locally uniformly in <C \ ([—1 — e, 1 + e]) for some e > 0. 

Proof. Obviously, the first statement is equivalent to the second one by construc- 
tion. So, the equivalence of (ii) and (iii) proves the theorem. First, suppose that 
(iii) holds true and, therefore, from Theorem 16.81 we see that (|6.1ip is fulfilled. 
Then, according to and (IT^ . the operator L is bounded. 

The implication from (ii) to (iii) follows from formulas (I4.2ip . (14.241) and Theo- 
rem |S3l □ 



In some cases, we can specialize Theorem [ 
Corollary 6.10. Let ^ admit the following representation 

j=i J 

where a is a finite nonnegative measure on [0, 1], tj < 0, Oj > for j = 1, . . . ,p, 
and S-i is a fixed real number. Then the following statements hold true: 

(i) the sequence {i?f"J^"^'}°^o °f diagonal Pade approximants converges to ^ 
locally uniformly in C\ ([0, 1] U {tk}^_^_i); 
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(ii) the polynomials Pj and Qj corresponding to a + J2^=i ~ ^j) possess 

the following property 



sup 



< oo; 



Qn,-l(0)-S-lPn,.-l(0) 

(iii) the operator £ in the factorization (j4.ip of J corresponding to 5" is bounded 

[0,! 



(iv) the operator £ — £(s_i) in the factorization (|4.15p of J with the parameter 
S-i = Sq corresponding to the function 

"1 da{t) " - 



t,- - A 



is bounded in l"^^ 



Proof. Clearly, ^ £ Dl^. Actually, statement (i) for the class in question was 
proved in 1101 . Statement (ii) is an immediate consequence of statement (i) and 
Theorem 16.81 The rest follows from Theorem 16.91 □ 



Remark 6.11. 1) In fact, in jTO] we proved the convergence results for the Pade 
approximants to Ns-functions having the asymptotic expansion (j2.1[) . 

2) As follows from (g^S]), (li^ the function F in (I5T^ belongs to the class S"", 
where k = p, ii F{0~) = /q^ ^ + 77 < 0, and k =p+l, if < F(O-) < 00. 

7. Examples 

7.1. The appearance of the block structure. In this subsection we give a 
family of Jacobi matrices such that the ChristofFel transform of these matrices 
leads to block tridiagonal matrices. 

Let us consider a sequence {PjyfLo of monic orthogonal polynomials satisfying 
the following three term recurrence relation 

AP,(A) = Pj+i(A) + Cj_iP,-_i(A), J e Z+, 

where Cj > and P_i(A) = 0, Po(A) = 1. For example, if cj = 1/4, j G Z+, 
then Pj are Chebyshev polynomials of second kind. Clearly, the polynomials Pj 
are symmetric, i.e. 

(7.1) P,-(-A) = (-l)^P,(A), jeZ+. 

Besides, the corresponding monic Jacobi matrix has the following form 

/O 1 \ 

Co 1 



J 



Cl 



and is associated with the Markov function 

dajt) 
_i i- A' 



P(A) 



where da is a symmetric measure. 

It follows from ([71]) that P2i+i(0) = and P2j(0) ^ for every j G Z+. The 
latter fact immediately implies that iXj = 2j, j G Z+, and, thus, ij — 2 for j e Z_|_. 



26 



MAXIM DEREVYAGIN AND VLADIMIR, DERKACH 



Now, Proposition [SUl gives the LC/-factorization J — LU, where L and U are block 
lower and upper triangular matrices having the forms 



Li It, 



L = 



\ 



L9 /fo 



V 



u = 



fUo Do 
Ui Di 

U2 



\ 



which entries are as follows 



L 



_ f C2J + 1 





1 

C2j 





1 



Multiplying L and U the other way around, we arrive at the monic generalized 
Jacobi matrix 

/So So \ 



UL = 



£1 ^2 



V 



where the entries have the following form 



1 

C2j 0^ 





Cj + lCj+2 



" 1 1 



Finally, observe that F can be represented as follows 

dp{t) 



F{\) = A 



e N. 



/o 

Thus, due to Theorem 13.21 we have that ,3 is associated with the function 



tdp{t) 



e d: 



7.2. Unboundedness of the Christoffel transform. Here we give an example 
of a bounded Jacobi matrix such that its Christoffel transform is an unbounded 
matrix. 

Consider the following monic 2-periodic Jacobi matrix (see [11] Example 1]) 



/ao 1 
1 ai 1 



1 02 



V 



associated with the function 
Fi\) 



A 1 / A(A2 - A-4) 
2 ^ 2V A^l ' 



where the branch is chosen so that F{X) as A 00 (see pT|). 
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Let US find the Lf/-factorization of J, that is, let us represent J as follows 

/ui 1 \ 
uih U2 + I1 



J = LU ■ 



1 



^2^2 ^3 + h 



\ 



I 



Clearly, we have the following relations 

u\ = 1, Ukh = 1, U2k + hk-i = 0, 'U2fe+i + /2fe = 1, fee N. 
Next, by induction we have that u\ = l,li = 1, U2 = —1, h = —1, and 

'"2fe+l = fc + l, U2/S+2 = --^7^, l2k+l = Y^, l2k+2 = -{k + 1), fceN. 

Therefore, the operators L and U in L?7-factorization of J are unbounded. More- 
over, the Christoffel transformation 



/wi +I1 1 
U2I1 U2 + h 



Jc = UL 



1 



\ / 2 1 

-1 -2 1 



W3^2 WS + ^3 



V 



-2 5/2 



V 



is also an unbounded monic generalized Jacobi matrix associated with the function 



.r.,.. . A2 , A /A{A2 - A-4) ^ 



[4; 
[s: 
[6: 
[7: 

[9 

[lo- 
in 

[12: 
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